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Abstract. We study the intermittent behavior of the energy decay and linear magnetic response of a glassy 
system during isothermal aging after a deep thermal quench using the Edward- Anderson spin glass model 
as a paradigmatic example. The large intermittent changes in the two observables are found to occur in a 
correlated fashion and through irreversible bursts, 'quakes', which punctuate reversible and equilibrium-like 
fluctuations of zero average. The temporal distribution of the quakes it foun to be a Poisson distribution 
with an average growing logarithmically on time, indicating that the quakes are triggered by record sized 
fluctuations. As the drift of an aging system is to a good approximation subordinated to the quakes, 
simple analytical expressions (Sibani et al. Phys Rev B 74, 224407, 2006) are available for the time and 
age dependence of the average response and average energy. These expressions are shown to capture the 
time dependencies of the EA simulation results. Finally, we argue that whenever the changes of the linear 
response function and of its conjugate autocorrelation function follow from the same intermittent events a 
fluctuation-dissipation-like relation can arise between the two in off-equilibrium aging. 

PACS. 65.60.+a Thermal properties of amorphous solids and glasses - 05. 40. -a Fluctuation phenomena, 
random processes, noise, and Brownian motion - 61.43.Fs Glasses - 75.10.Nr Spin-glass and other random 
models 



1 Motivation 

In noise spectra from mesoscopic aging systems, reversible 
fluctuations are punctuated by rare and large, so called 
intermittent, events [?,?,?,?,?,?], which arguably signal 
switches from one metastable configuration to another [?,?]. 
The intermittent events usually appear as an exponen- 
tial tail in the Probability Density Function (PDF) of 
the fluctuations. The central part of the PDF describes 
equilibrium-like behavior with its zero-centered Gaussian 
shape. Exploiting this information can adds new twists 
to long debated issues as the multiscale nature of glassy 
dynamics and the associated memory behavior. 

A record dynamics scenario for aging [?,?,?] builds on 
two main assumptions: (i) record sized energy fluctuations 
within metastable domains trigger irreversible intermit- 
tent events, or quakes; (ii) these, in turn, control all sig- 
nificant physical changes, e.g. heat release, configurational 
decorrelation and linear magnetic response. In brief, all 
important changes arc considered to be subordinated to 
the quakes, with the latter triggered by energy fluctua- 
tions of record size. Combining the information from the 
equilibrium-like fluctuations with non-thermal properties, 
e.g. irreversible energy losses, leads to testable predictions 
for the time and temperature dependencies of the fluctu- 
ation spectra which have been confirmed in a number of 
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While macroscopic (average) linear response functions 
have long been the tool of choice for probing aging in 
magnetic systems, [?,?,?,?,?], the predictions of record 
dynamics have so far only been tested for experimental 
Thermoremanent Magnetization (TRM) data, spin glass 
data [?] . Simulations offer certain advantages over experi- 
ments: it is possible to simulate an instantaneous quench, 
simply by choosing a 'random', i.e. high temperature, ini- 
tial condition. Notably, an instantaneous initial quench 
leads to the 'full aging' [?], scaling behavior simply ex- 
plained in record dynamics; secondly, the statistical anal- 
ysis is considerably simplified when thousands of indepen- 
dent traces can be generated; lastly, temporal correlations 
between the intermittent changes of the magnetization 
and the energy can be extracted. 

The Edwards- Anderson (EA) spin-glass model used in 
this work is a paradigmatic example of an aging system. 
Its macroscopic (average) magnetic response and autocor- 
relation decay have been thoroughly investigated [?,?,?,?] 
and its mesoscopic fluctuations properties have attracted 
some recent attention [7,7,7], Here, data from extensive 
simulations are analyzed with focus on the intermittency 
of the energy and Zero Field Cooled Magnetization (ZFCM) 
fluctuations We confirm that quakes carry the net drift 
of the energy [7] and correlate strongly with the large 
intermittent magnetization fluctuations carrying the net 
change of the linear response. The idea [7, 7] that quakes 
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Fig. 1. (Color on line) The left panel shows the PDF (blue circles) of the magnetization changes 8M occurring over intervals 
of length St = 20 during isothermal simulations at temperature T = 0.4. The statistics is based on 2 • 10 5 data points obtained 
in the age interval [1000, 3000] using 2000 independent trajectories. The zero average Gaussian (black line) fits data points with 
—80 < 5M < 50. The right panel shows the PDF of the energy changes (blue circles), and a Gaussian fit of the positive energy 
fluctuations (black line). The red squares pertain to the subset of energy changes occurring in the immediate vicinity of large 
[5M > 50 ) magnetization changes. All data are obtained from the same set of trajectories as those shown in the left panel. 



have a Poisson distribution with a logarithmic time depen- 
dence is derived from record statistics and is central to the 
theory. To verify it empirically, we consider the temporal 
statistics of the difference of 'logarithmic waiting times' 
log(tg) — log(f 9 _i), where t q marks the occurence of the 
g'th quake in a given trace. For the above Poisson distri- 
bution, these logarithmic differences are exponentially dis- 
tributed. General mathematical arguments lead to eigen- 
value expansions for the dependence of the average energy 
and linear response on the number of quakes, and then, 
via the subordination hypothesis, to power-law expansions 
for the time dependences of the same quantities [?,?]. The 
expansions are tested below against the EA model nu- 
merical simulation results. The origin of approximate off- 
equilibrium Fluctuation-Dissipation like relations is dis- 
cussed in the last section from the point of view of record 
dynamics. 

Finally, a notational issue: as in ref. [?], the variable 
t generically denotes the time elapsed from the initial 
quench, i.e. the system age. The external field is switched 
on at time t = t w , and t b s — t — t w denotes the obser- 
vation time, a quantity called t in refs. [?,?] and in many 
experimental papers e.g. [?,?]. Unless otherwise stated, 
we denote the average energy and magnetization by [i e 
and fizFCM) and reserve the symbols e and M for the 
corresponding fluctuating quantities measured in the sim- 
ulations. The simulation temperature is denoted by T. 



2 Model and Simulation methods 

In the Edwards- Anderson model, N Ising variables, <x; = 
±1, are placed on a cubic lattice with toroidal boundary 
conditions. Their interaction energy is 

N N 

E = -^^o-idjJij - Hj2°i, (1) 

i,j i 

where J^- <7j = M is the magnetization and H is the mag- 
netic field. The non-zero elements of the symmetric inter- 
action matrix Jij connect neighboring sites on the lattice. 
The interactions are drawn from a Gaussian distribution 
with zero average and unit variance, a choice setting the 
scale for both temperature and magnetic field. 

In the simulations, we use system size N = 16 3 and 
collect the statistics of energy and magnetization changes 
from either several thousands of independent trajectories, 
each corresponding to a different realization of the J^'s. 
After the initial quench, aging procedes isothermally in 
zero field untill time t — t w , where a small magnetic 
field H = 0.1 is instantaneously turned on. Isothermal 
simulations lasting up to age t = 100 t w are carried out 
for t w = 50, 100, 200, 500 and 1000 at temperatures T = 
0.1,0.2, . . .0.8. The simulation engine utilizes the Waiting 
Time Method [?,?], a rejectionless or 'event driven' algo- 
rithm which operates with an 'intrinsic' time, loosely cor- 
responding to a sweep of the Metropolis algorithm. Unlike 
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Fig. 2. (Color on line) The left panel depicts the PDF of the magnetization changes SM over intervals of length St. The 
simulation temperature is T = 0.4. Green circles, blue squares and magenta diamonds correspond to St — 20, 40 and 60, 
respectively. The data are collected in the interval [1000, 3000]. In order of increasing St, the statistics is based on 10 5 , .5 -10 s 
and .33 • 10 data points. The full line is a fit of all data with —80 < 5M < 50 to a Gaussian centered at zero. For 5M > 50 the 
PDF clearly deviates from a Gaussian shape. The intermittent tail grows with St and fully determines the average change in 
the ZFC magnetization. The insert shows dM, averaged over t b 3 , versus St. The right panel shows the variance of the (fitted) 
Gaussians as a function of the reduced temperature T/T g . 



the Metropolis algorithm, the WTM can follow spatio- 
temporal patterns on very short time scales, a property 
advantageous for intermittency studies. 

3 Simulation results 

In this Section, simulation results are presented together 
with some relevant theoretical considerations. The first 
subsection deals with the distribution of energy and mag- 
netization changes, SE and SM, occurring over small inter- 
vals St. The Probability Density Function (PDF) of these 
quantities characterizes intermittency in the EA model in 
a simple and direct way. In the second subsection, the 
quakes are (approximately) indentified within each data 
stream, and the temporal aspects of their distribution are 
studied. In the third and last subsection, the average lin- 
ear response and average energy obtained from the simu- 
lations are compared to analytical predictions which are 
based on the temporal distribution of the quakes. 

3.1 Energy and magnetization fluctuation statistics 

For a first statistical description of the energy and the 
ZFCM fluctuations, energy and magnetization changes SE 
and SM occurring over small intervals St are sampled dur- 
ing the time interval [1000, 3000]. The left panel of Fig. [1] 



shows the PDF of the magnetic fluctuations (blue cir- 
cles) on a log scale. The (black) line is the fit to a zero- 
centered Gaussian obtained using data within the interval 
—80 < SM < 50. The distribution of the remaining large 
positive magnetization changes is seen to have exponential 
character. The same runs are used to collect the PDF of 
the energy fluctuations, shown in the right panel of Fig. [1] 
(lower data set, blue circles). Again, we see a combination 
of a zero centered Gaussian and an intermittent exponen- 
tial tail. The full line is a fit to the zero centered Gaussian 
in the interval -5 < SE < 20. 

As intermittent fluctuations are rare, the temporal cor- 
relation between intermittent energy and ZFC magnetiza- 
tion changes is best observed using conditional PDF's (up- 
per curve, red squares). The conditional PDF shown only 
includes those SE values which either fall in the same or 
in the preceding small interval of width St as magnetic 
fluctuations above the threshold SM = 50. In the left 
panel of the figure, the threshold is seen to be near the 
boundary between reversible and intermittent magnetiza- 
tion changes. Note how the Gaussian parts of the full and 
of the conditional PDFs nearly coincide, while the cor- 
responding intermittent tails differ. The tail of the condi- 
tional PDF is strongly enhanced (roughly, between 10 and 
100 times) when only fluctuations near a large magnetiza- 
tion change are included in the count. A similar situation 
is observed at any low temperatures. 
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Fig. 3. (Color on line) For eight different temperatures T — 
0.2, 0.3 . . . 0.9 (circles, dots, squares, plusses, diamonds, stars, 
pentagrams, and hexagrams) the estimated cumulative prob- 
ability Prob(log tq/tq-i > x) is plotted versus x, with t q , 
q — 1,2... denoting the times at which the quakes occur. The 
quakes are identified as described in the main text using the 
filtering parameter value / = 7. The dotted lines are fits to 
the exponential exp(— ax). The collapse of the eight data sets 
demonstrates that the temperature dependence is very weak, 
except at the lowest T. In the insert, the fitted a values are 
plotted versus T. The circles pertain to the data shown in the 
main plot, and the squares pertain to corresponding data ob- 
tained using / = 8 as a filter parameter. 



Figure [2] details the dependence of the PDF of the 
magnetization fluctuations on St and T: The three PDFs 
shown in the left panel of Fig. [2] as green circles, blue 
squares and magenta diamonds are obtained in the inter- 
val [1000,3000], for T = 0.4 and for St = 20,40 and 60, 
respectively. The black line is a zero centered Gaussian fit- 
ted to the PDF in the range -80 < SM < 10. The fit is ob- 
tained by optimizing the variance of the Gaussian distribu- 
tion, Oq. The latter quantity should not be confused with 
the variance of the full distribution, which is much larger 
and heavily influenced by the tail events excluded from 
the Gaussian fit. Importantly, since the Gaussian part of 
the PDF is independent of St quasi-equilibrium ZFCM 
fluctuations at T = 0.4 can be treated as uncorrelated for 
times larger than St = 20. Secondly, since the Gaussian is 
centered at zero, any changes in the average magnetiza- 
tion are exclusively due to the intermittent events. Similar 
conclusions were reached for the spin-glass TRM magneti- 
zation, [?] and for the energy outflow in the EA model, [?] 
and in a p-spin model with no quenched randomness. [?] 
The insert describes the dependence on St of the magneti- 
zation change (SM) averaged over the observation interval 
[1000, 3000]. E.g. the data point corresponding to St = 20 
depicts the average of the lowest PDF plotted in the main 
figure. The right panel of fig. [2] shows a plot of the ratio 



Oq/T versus T for the Gaussian part of the magnetic fluc- 
tuations. As the Fluctuation-Dissipation Theorem (FDT) 
applies to the equilibrium-like part of the dynamics, the 
ordinate can be interpreted as the (gedanken) linear mag- 
netic susceptibility of metastable configurations. The over- 
all shape of the curve is reminiscent of the experimental 
T dependence of the ZFCM below T g . [?] 



3.2 Temporal distribution of intermittent events 

The PDF of intermittent fluctuations in the EA model 
scales with the ratio St/t, where t is the age of the sys- 
tem at the beginning of the sampling interval [?]. The 
rate of quakes accordingly decays as the inverse of the 
age, in agreement with the claim that the the number 
of independent quakes in an interval (t, t') has a Pois- 
son distribution with average (nj) = alog(t/t') The pa- 
rameter a characterizing the average is interpreted as the 
number of thermalized domains contributing in parallel 
to the fluctuation statistics. This quantity is expected to 
be temperature independent and linearly dependent on 
system size [?]. Using that the differences r q = log(i g ) — 
log(t 9 _i) = log(tq/t q -x) are independent random num- 
bers, all exponentially distributed with the same average 
(Tq) = 1/a, identifying the times t q , q = 1, 2, . . . at which 
the quakes occur allows one to check the temporal aspects 
of statistics (as opposed to the distributionof quake sizes) . 

The identification entails some challanges of statisti- 
cal and patter-recognition nature. The residence times 
t q —t q -\ expectedly have a broad distribution [?], which, in 
connection with a finite sampling time creates a negative 
sampling bias on large t q values [?] , even with the quakes 
unabiguously identified. Secondly, correlations may well 
be present between large and closely spaced spikes in the 
signal which collectively represent a change of attractor. 
Last but not least, identification of the t q 's within a trace 
must rely on the negative sign of the fluctuations and 
on their 'sufficiently large' size. I.e. by definition suffi- 
cient energy must be released to make the reverse process 
highly unprobable within an observation span stretching 
up to 2t for a quake which occurs at t. [?] A related and 
generic property of aging systems (also directly accessi- 
ble by an intermittency analysis [?]) is that reversible 
equilibrium-like fluctuations dominate the dynamics on 
time scales shorter than the age t. A simple empirical cri- 
terion for distinguishing reversible fluctuations from irre- 
versible quakes uses an upper bound on the probability, 
P rev (t, SE) that at least one thermal energy fluctuation 
of positive sign be among the t/St observations. For small 
exp(—SE/T), this happens with probability P rev {t, SE) ss 
t/Stcexp(—SE/T), where the positive number c is un- 
known. A negative energy fluctuation is labeled as a quake 
if its absolute value \SE\ satisfies P rev (t,\SE\) < 10 _n , 
where n is a positive number. For a thermal fluctuation 
which occurs at age t, this criterion produces the inequal- 
ity 
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>/ + logt; where / = nln(10) + ln(c/<ft). 

(2) 
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The formal definition of the 'filter' parameter / on the 
right hand side of the inequality contains the unknown 
parameter c as well as a free parameter n. We must there- 
fore treat / as a free parameter regulating how strict a 
filter the fluctuations must pass to qualify as quakes. 

For temperatures T — 0.2, 0.3 . . . 0.9, a range spanning 
most of the EA spin-glass phase, simulations were per- 
formed in the interval t S (200, 10200). Within each tra- 
jectory, quakes are identified by applying Eq. [2j and the 
quantities r q = log(t q /t q -i) are calculated and binned. 
The procedure is repeated for 2000 independent trajec- 
tories for each value of the temperature. The magnetic 
field values H — 0.1 for t > t w = 200 and zero other- 
wise, the same as in all other simulations, were used for 
simplicity. Note however that, as also demonstrated by 
Fig. 21 the magnetic field has no discernible effects on 
the energy statistics. The standard deviation, er(i), of the 
empirical probability of a point falling in the i'th bin is 
cr ( i ) = \f{Pth{i)/N obs ), where p t h(i) is the corresponding 
(unknown) theoretical probability. N b s is the number of 
r q values collected. Replacing p t h with the corresponding 
empirical probability, yields an estimate for a(i). The em- 
pirical cumulative distribution Prob(r g = log(t q /t q -i) > 
x) is fitted to the exponential be ax , with fit parameters 
obtained by mimizing the sum of the square differences 
to the data points, each weighted by the reciprocal of the 
(estimated) variance of the data point. The prefactor b has 
no physical significance. It would be unity, where not for 
the fact that the lowest x value, for which the ordinate is 
equal to one, is located at half bin size, rather than at zero. 
To check the influence of / on the form of the distribu- 
tion, all the simulations were done twice, using the values 
/ = 7 and / = 8. The results shown in the main panel of 
Fig. [3] pertain to the / = 7 case. The results for / = 8 are 
very similar, i.e. both cases produce an exponential dis- 
tribution. The scale parameter a lacks any significant T 
dependence, except for the lowest temperature, and seems 
only weakly / dependent. From the statistics, the number 
of domains can be estimated as a « 20. Dividing the total 
number of spins by a an order of magnitude estimate is 
obtained for the number of spins in a thermalized domain. 
The corresponding linear size comes out around 6 spins, a 
figures which compares well to the range of domain sizes 
observed in domain growth studies of the EA model [?]. 
The scale invariance of the energy landscape implict in 
the T independence of a is to a large degree confirmed. 
It is expected that the largest deviations be found at the 
lowest temperatures, where the discreteness of the energy 
spectrum begins to make itself felt. 

3.3 Average energy and linear response 

The drift of the average energy /i e and average magne- 
tization fizFCM was seen to mainly depend on the the 
number of quakes nj which fall in the relevant observa- 
tion interval. For an interval (t, t') this number was shown 
to have a Poisson distribution with average 

(n/) = a\n(t/t'). (3) 



For a final check of the applicability of record dynamics we 
now compare available analytic formulas [?,?] for the av- 
erage energy, magnetic response and magnetic correlation 
function to the spin glass data. The subordination hypoth- 
esis used to derive the formulas fully neglects the effects 
of pseudo-equilibrium fluctuations. Then, all observables 
have a rather simple and generic nj dependence — a su- 
perposition of exponential functions. Averaging over the 
distribution of nj produces a superposition of power-laws 
terms, some of which can be further expanded into a log- 
arithmic and a constant term. The independent variables 
are t and t/t w for the energy and the magnetization, re- 
spectively, since all quakes contribute to the energy decay, 
while only those falling between t w and t contribute to the 
response. 

Because pseudo-equilibrium fluctuations are excluded, 
deviations between predicted and observed behavior are 
to be expected. Indeed, as discussed later, the linear re- 
sponse, data have a small additive t w dependence which 
is beyond the reach of the description and which is in 
most studies explicitely split off as a 'stationary contribu- 
tion [?]. 

The average energy and linear response are calculated 
using 4000 independent trajectories. The average energy 
/i e (i) is found to decay toward an (apparent) asymptotic 
limit ec, according to the power-law 

Me(*) - £G = a e t K . (4) 

In the first five panels of Fig. [4j the full line corresponds to 
Eq.[¥]and the symbols are empirical estimates of the aver- 
age energy. Each panel shows data obtained at the temper- 
ature indicated. For each T, the five data sets displayed are 
for t w = 50, 100, 200, 500 and 1000 (green, blue, magenta, 
black and red, respectively). The lack of a visible t w depen- 
dence shows that the magnetic contribution to the energy 
is negligible, as expected. The last panel of Fig. [4] summa- 
rizes the dependence of ec and of the decay exponent A e on 
the reduced temperature T/T c . The estimate T c — 0.98 [?] 
is used for the critical temperature of the EA model. In- 
terestingly, the fitted value of eq is nearly independent of 
T, and is close to the estimated ground state energy of 
the EA model, [?] e G = -1.7003 ± 0.008. Eq. H concurs 
with the observation that a 'putative' ground state en- 
ergy of complex optimization problems can be guessed at 
early times. [?,?] The prefactor a e has a negligible T de- 
pendence. In contrast, the exponent A e has a clear linear 
dependence A e = —0.05 — 0.38^. 

For \ e \nt < 1, the average energy decay is logarith- 
mic, He(t) ea + a e + a e X e ln(t) . . ., and the decay rate 
falls off as the inverse of the age [?,?]. The coefficient a e X e , 
which receives its temperature dependence from A e , is pro- 
portional to the average size of a quake [?] . Even though 
quakes are strongly exothermic, their average size may 
slightly increase with T. This effect of thermal activation 
on the size of quakes is analyzed (for a different model) in 
Ref. [?]. 
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Fig. 4. (Color on line) The first five panels depict, for the temperatures indicated, scaling plots of the energy per spin, with 
the constant ec, subtracted. The values of t range from i = t w to 400t w , with t w — 50,100,200,500 and 1000, (green, blue, 
magenta, black, and red circles respectively). The black line is a fit to EqU The value of t w is irrelevant for the energy decay. 
The last panel shows the T dependence of the parameters of Eq.|4]on the reduced temperature T/T c . 
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Fig. 5. (Color on line) The first five panels display, for the temperatures indicated, scaling plots of the magnetization with a 

small t w dependent quantity b(t w ) added. The abscissa is t/t w , with t w equal to 50, f00, 200, 500 and 1000. The values of t range 

from t — t w to t — 400^. The data are fitted to the sum of a logarithm and a power-law (see Eq(5}. In the inserts, the relaxation 

rate S is plotted versus the observation time t — t w . The values of S are obtained by applying Eq. [6] to the fitted functional 

form for \izfcm given in Eq[5] The last panel shows the dependence of the parameters of Eq. [5] on the reduced temperature 
t It tu„ + „f \ ;„ „«.„«u;„„l „„„„„„„ 



8 Paolo Sibani: Linear response in aging glassy systems, intermittency and the Poisson statistics of record fluctuations. 



An excellent parameterization of the average ZFCM 
time dependence is given by 

HZFCAt(t/t w ) =b + a m In (j~^J + b m ' ( 5 ) 

In the theory, the left-hand side should only depend on 
t/t w , i.e. b should be constant. However, a small t w de- 
pendent term, the well known stationary part of the re- 
sponse, is present. To be able to nevertheless plot the data 
versus t/t w we shift them vertically in a t w dependent fash- 
ion, i.e. we subtract the stationary term. The magnitude 
of the shift increases with t w and with T, reaching, at its 
highest, approximately 15% of the value of 60 which fits 
the tw — 50 data. 

The first five panels of Fig. [5] show the average ZFCM 
(dots) and the corresponding fits (lines) as a function of 
t/t w , for t w = 50,100,200,500 and 1000. The simula- 
tion temperature used is indicated in each case. Sets of 
data corresponding to different t w values are color coded 
as done for the average energy. The black line is given 
by Eq. with its three parameters determined by least 
square fits. The field switch at t w has a 'transient' ef- 
fect described by the power-law decay term: When t/t w 
is sufficiently large, the ZFCM increases proportionally to 
ln(t/t w ), and hence at a rate decreasing as a m /t. In the 
last panel of Fig. the decay exponent A m and the con- 
stants a m and b m are plotted versus T. The full line shows 
the linear fit — A m = 1.10 + 0.33T/T C , the dotted lines are 
guides to the eye. Importantly, the constant of propor- 
tionality a m is practically independent of T, a property 
also shared by a e . This strongly indicates that the loga- 
rithmic rate of quakes, a (see Eq. [3]) must similarly be T 
independent. 

The derivative of hzfcm with respect to the logarithm 
of the observation time, the so-called 'relaxation rate', can 
be written as 

S(t w ,t) = (t - t w ) d ^ Z ^ CM =(t- t w )r Z FCM, (6) 
at 

where Tzfcm l& the rate of magnetization increase. The 
plots of S versus the observation time t b s — t t^ are 
produced via Eqs. [5] and and dispayed in the inserts of 
the first five panels of Fig. [5l One recognizes the char- 
acteristic peak at t b s = t w [?,?]. The limiting value for 
t 3> t w is a m , precisely the asymptotic logarithmic rate of 
increase of the ZFC magnetization.. 

The sum of the TRM and ZFCM is the field cooled 
magnetization (FCM). Unlike the ZFCM, both FCM and 
TRM have a a non-linear H dependence. [?] The latter 
does not seem to affect the intermittent behavior, at least 
if numerical simulation data and experimental data can 
be treated on the same footing: Experimental TRM de- 
cay data [?] have the same general behavior as the ZFCM 
data just discussed. However two power-law terms, rather 
than one, are needed to describe the, more copmlex, ex- 
perimental transient. The exponents have ranges similar 
to A m , and a somewhat stronger variation with T. The 
large t/t w asymptotic behavior is a logarithmic decay, 



aln(t/t w ), where, importantly, the coefficient of propor- 
tionality a is T independent, except very close to T g . 

With an eye to the following discussion, we recall that 
the average configuration autocorrelation function of the 
EA model, normalized to unity at t = t w , decays between 
times t w and t as 

Hc(t/t w ) = {t/t w ) Xc ; t > t w , (7) 

where the exponent can be fitted to the linear T depen- 
dence X C (T) = —0.25^. [?] The algebraic decay follows 
from the record dyanmics. As for the response, the the- 
oretical arguments neglect the effect of quasi-equilibrium 
fluctuations for t « t w . Furthermore, they fail near the fi- 
nal equilibration stages. Note that in the notation of Ref. [?] 
the ratio t/t w is written as 1 + t/t w . 

4 Discussion and Conclusions 

In this and a series of preceding papers [?,?,?,?,?,?], we 
have argued that non-equilibrium events, the quakes, are 
key elements of intermittency and non-equilibrium ag- 
ing. The approach takes irreversibility into account at 
the microscopic level, stressing that thermal properties 
alone, e.g. free energies, are inadequate to explain non- 
equilibrium aging. These properties remain nevertheless 
important in the description, since the record sized (posi- 
tive) energy fluctuations which are assumed to trigger the 
quakes are drawn from an equilibrium distribution. As 
earlier discussed, [?] the energy landscape of each domain 
must be self-similar in order to support a scale invariant 
statistics of record fluctuations. In this optics, the scale in- 
variance of the 'local' energy landscape attached to each 
domain, rather than the scale invariance of real space ex- 
citations, is at the root the slow relaxation behavior of 
aging systems [?,?,?]. 

A different approach to non-equilibrium aging consid- 
ers a generalization of the Fluctuation Dissipation The- 
orem (FDT) [?,?,?,?]. Its relation to record-dynamics is 
briefly considered in below. Out of equilibrium, the FDT 
is never fulfilled exactly [?] , nor is it generally possible to 
write the linear response as a function of the conjugate au- 
tocorrelation alone. Nevertheless, for time scales t — t w << 
t w , the drift part of aging dynamics is negligible, and the 
FDT does, in practice, apply [?]. As mentioned, conju- 
gate response and autocorrelation functions are naturally 
divided into stationary and non- stationary parts, which 
pertain to the pseudo-equilibrium and off-equilibrium ag- 
ing regimes, respectively. Adopting the standard notation 
C for the correlation and R for the magnetic response, 
the FDT reads (1 - C(t))/T = R(t)/H. For equilibrium 
data, a plot of C versus x = R/H yields a straight line 
with slope — 1/T. For aging data, the same plot produces 
a straight line in the quasi-equilibrium regime, e.g. at 
early times where C is close to one. Accordingly, a mea- 
sure of the deviation from quasi-equilibrium is the Fluc- 
tuation Dissipation ratio (FDR) [?], which is defined as 
X(C) = -Tdx(C)/dC, In the case of the EA model, the 
relation InC = A c log(t/t w ), which follows from inverting 
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Eq. [JJ can be used to find x(C) an d to calculate X(C). 
The latter is nowhere a linear function of C, which is ex- 
pected, as the stationary fluctuation regime is excluded 
from the description. The effective temperature [?] is usu- 
ally defined from the large t and t w asymptotic limit of the 
Fluctuation Dissipation ratio. Effective temperatures may 
depend on the choice of conjugate observables [?], and are 
not easily measured experimentally [?,?], but offer never- 
theless a simple and appealing characterization of aging 
dynamics. 

In a record-dynamics context, fluctuation-dissipation 
like relations arise out of equilibrium because correlation 
and response are both subordinated to the same quakes. 
Due to the monotonicity of their t/t w dependencies, each 
of the two can can be written as a function of the other, 
and a FDR can thus be constructed. Asymptotically, both 
correlation and response may have an approximate log- 
arithmic dependence on t/t w : From Eq. [5l we see that 
R{t/t w ) a m ln(t/t w ) for t/t w 3> 1. Considering that 
A c vanishes for T — > 0, we also obtain, for the same 
range of t/t w values, the inequality |A C (T)| hx(t/t w ) < 1. 
When the inequality is fulfilled, Eq. [7J can be written as 
1 - C(t/t w ) « -A C (T) \n(t/t w ) cx R(t/t w ). A relation for- 
mally similar to the FDT holds, with T replaced by 

_ -\ C H _ T 
T »tt — — 77777 7777 ■ ( 8 ) 

As a m (x H, there is of course no H dependence. Im- 
portantly, T c ff oc T, clS Oifri IS independent of T. For the 
parameter values obtained from the fits, T e ff > T, even 
though a general argument to support the inequality is 
lacking at the moment. 

The approach used in this paper and in ref. [?] should 
be generally applicable to check for the presence of record- 
dynamics features in intermittent fluctuation data. Clearly 
the ability to perform calorimetry experiments is essential 
to directly check the temperal statistics of the quakes. In 
the absence of such data, one may assume that quakes 
produce large intermittent changes in other observables, 
and use their statistics instead. 
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Abstract. We study the intermittent behavior of the energy decay and linear magnetic response of a glassy 
system during isothermal aging after a deep thermal quench using the Edward- Anderson spin glass model 
as a paradigmatic example. The large intermittent changes in the two observables are found to occur in a 
correlated fashion and through irreversible bursts, 'quakes', which punctuate reversible and equilibrium-like 
fluctuations of zero average. The temporal distribution of the quakes it foun to be a Poisson distribution 
with an average growing logarithmically on time, indicating that the quakes are triggered by record sized 
fluctuations. As the drift of an aging system is to a good approximation subordinated to the quakes, 
simple analytical expressions (Sibani et al. Phys Rev B 74, 224407, 2006) are available for the time and 
age dependence of the average response and average energy. These expressions are shown to capture the 
time dependencies of the EA simulation results. Finally, we argue that whenever the changes of the linear 
response function and of its conjugate autocorrelation function follow from the same intermittent events a 
fluctuation-dissipation-like relation can arise between the two in off-equilibrium aging. 

PACS. 65.60.+a Thermal properties of amorphous solids and glasses - 05. 40. -a Fluctuation phenomena, 
random processes, noise, and Brownian motion - 61.43.Fs Glasses - 75.10.Nr Spin-glass and other random 
models 



1 Motivation 

In noise spectra from mesoscopic aging systems, reversible 
fluctuations are punctuated by rare and large, so called 
intermittent, events |1I2I3I4I5[q] . which arguably signal 
switches from one metastable configuration to another [7|8j . 
The intermittent events usually appear as an exponen- 
tial tail in the Probability Density Function (PDF) of 
the fluctuations. The central part of the PDF describes 
equilibrium-like behavior with its zero-centered Gaussian 
shape. Exploiting this information can adds new twists 
to long debated issues as the multiscale nature of glassy 
dynamics and the associated memory behavior. 

A record dynamics scenario for aging [8 9 10J builds 
on two main assumptions: (i) record sized energy fluctua- 
tions within metastable domains trigger irreversible inter- 
mittent events, or quakes; (ii) these, in turn, control all 
significant physical changes, e.g. heat release, configura- 
tional decorrelation and linear magnetic response. In brief, 
all important changes are considered to be subordinated 
to the quakes, with the latter triggered by energy fluctu- 
ations of record size. Combining the information from the 
equilibrium-like fluctuations with non-thermal properties, 
e.g. irreversible energy losses, leads to testable predictions 
for the time and temperature dependencies of the fluctu- 
ation spectra which have been confirmed in a number of 
cases |8lllll2ll3IT4ll5llri| . 



While macroscopic (average) linear response functions 
have long been the tool of choice for probing aging in mag- 
netic systems, |17ll8ll9l2QI21j . the predictions of record 
dynamics have so far only been tested for experimental 
Thermoremanent Magnetization (TRM) data, spin glass 
data [16) . Simulations offer certain advantages over exper- 
iments: it is possible to simulate an instantaneous quench, 
simply by choosing a 'random', i.e. high temperature, ini- 
tial condition. Notably, an instantaneous initial quench 
leads to the 'full aging' [20], scaling behavior simply ex- 
plained in record dynamics; secondly, the statistical anal- 
ysis is considerably simplified when thousands of indepen- 
dent traces can be generated; lastly, temporal correlations 
between the intermittent changes of the magnetization 
and the energy can be extracted. 

The Edwards- Anderson (EA) spin-glass model used in 
this work is a paradigmatic example of an aging system. 
Its macroscopic (average) magnetic response and autocor- 
relation decay have been thoroughly investigated [22 2 3|24|25j 
and its mesoscopic fluctuations properties have attracted 
some recent attention |26l8ll3j . Here, data from exten- 
sive simulations are analyzed with focus on the intermit- 
tency of the energy and Zero Field Cooled Magnetiza- 
tion (ZFCM) fluctuations We confirm that quakes carry 
the net drift of the energy [8] and correlate strongly with 
the large intermittent magnetization fluctuations carry- 
ing the net change of the linear response. The idea [10|8j 
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Fig. 1. (Color on line) The left panel shows the PDF (blue circles) of the magnetization changes 8M occurring over intervals 
of length St = 20 during isothermal simulations at temperature T = 0.4. The statistics is based on 2 • 10 5 data points obtained 
in the age interval [1000, 3000] using 2000 independent trajectories. The zero average Gaussian (black line) fits data points with 
—80 < 5M < 50. The right panel shows the PDF of the energy changes (blue circles), and a Gaussian fit of the positive energy 
fluctuations (black line). The red squares pertain to the subset of energy changes occurring in the immediate vicinity of large 
[5M > 50 ) magnetization changes. All data are obtained from the same set of trajectories as those shown in the left panel. 



that quakes have a Poisson distribution with a logarith- 
mic time dependence is derived from record statistics and 
is central to the theory. To verify it empirically, we con- 
sider the temporal statistics of the difference of 'logarith- 
mic waiting times' \og{t q ) — log(t ? _i), where t q marks the 
occurence of the g'th quake in a given trace. For the above 
Poisson distribution, these logarithmic differences are ex- 
ponentially distributed. General mathematical arguments 
lead to eigenvalue expansions for the dependence of the av- 
erage energy and linear response on the number of quakes, 
and then, via the subordination hypothesis, to power-law 
expansions for the time dependences of the same quanti- 
ties [13I16J . The expansions are tested below against the 
EA model numerical simulation results. The origin of ap- 
proximate off-equilibrium Fluctuation-Dissipation like re- 
lations is discussed in the last section from the point of 
view of record dynamics. 

Finally, a notational issue: as in ref. [16j . the vari- 
able t generically denotes the time elapsed from the initial 
quench, i.e. the system age. The external field is switched 
on at time t = t w , and t t s = t — t w denotes the ob- 
servation time, a quantity called t in refs. |8|13j and in 
many experimental papers e.g. [27120) . Unless otherwise 
stated, we denote the average energy and magnetization 
by \x e and [Izfcm, and reserve the symbols e and M for 
the corresponding fluctuating quantities measured in the 
simulations. The simulation temperature is denoted by T . 



2 Model and Simulation methods 

In the Edwards- Anderson model, N Ising variables, crj 
■ I . are placed on a cubic lattice with toroidal boundary 
conditions. Their interaction energy is 

N N 

E=--Y^aiPjJij-H^2,a u (1) 

i,j i 

where J^- <7j = M is the magnetization and H is the mag- 
netic field. The non-zero elements of the symmetric inter- 
action matrix Jij connect neighboring sites on the lattice. 
The interactions are drawn from a Gaussian distribution 
with zero average and unit variance, a choice setting the 
scale for both temperature and magnetic field. 

In the simulations, we use system size N = 16 3 and 
collect the statistics of energy and magnetization changes 
from either several thousands of independent trajectories, 
each corresponding to a different realization of the J^'s. 
After the initial quench, aging procedes isothermally in 
zero field untill time t — t w , where a small magnetic 
field H = 0.1 is instantaneously turned on. Isothermal 
simulations lasting up to age t = 100 t w are carried out 
for t w = 50, 100, 200, 500 and 1000 at temperatures T = 
0.1, 0.2, ... 0.8. The simulation engine utilizes the Wait- 
ing Time Method |28|15j . a rejectionless or 'event driven' 
algorithm which operates with an 'intrinsic' time, loosely 
corresponding to a sweep of the Metropolis algorithm. Un- 
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Fig. 2. (Color on line) The left panel depicts the PDF of the magnetization changes SM over intervals of length St. The 
simulation temperature is T = 0.4. Green circles, blue squares and magenta diamonds correspond to St — 20, 40 and 60, 
respectively. The data are collected in the interval [1000, 3000]. In order of increasing St, the statistics is based on 10 5 , .5 -10 s 
and .33 • 10 data points. The full line is a fit of all data with —80 < 5M < 50 to a Gaussian centered at zero. For 5M > 50 the 
PDF clearly deviates from a Gaussian shape. The intermittent tail grows with St and fully determines the average change in 
the ZFC magnetization. The insert shows dM, averaged over t b 3 , versus St. The right panel shows the variance of the (fitted) 
Gaussians as a function of the reduced temperature T/T g . 



like the Metropolis algorithm, the WTM can follow spatio- 
temporal patterns on very short time scales, a property 
advantageous for intermittency studies. 

3 Simulation results 

In this Section, simulation results are presented together 
with some relevant theoretical considerations. The first 
subsection deals with the distribution of energy and mag- 
netization changes, SE and SM, occurring over small inter- 
vals St. The Probability Density Function (PDF) of these 
quantities characterizes intermittency in the EA model in 
a simple and direct way. In the second subsection, the 
quakes are (approximately) indentified within each data 
stream, and the temporal aspects of their distribution are 
studied. In the third and last subsection, the average lin- 
ear response and average energy obtained from the simu- 
lations are compared to analytical predictions which are 
based on the temporal distribution of the quakes. 

3.1 Energy and magnetization fluctuation statistics 

For a first statistical description of the energy and the 
ZFCM fluctuations, energy and magnetization changes SE 
and SM occurring over small intervals St are sampled dur- 
ing the time interval [1000, 3000]. The left panel of Fig. [1] 



shows the PDF of the magnetic fluctuations (blue cir- 
cles) on a log scale. The (black) line is the fit to a zero- 
centered Gaussian obtained using data within the interval 
—80 < SM < 50. The distribution of the remaining large 
positive magnetization changes is seen to have exponential 
character. The same runs are used to collect the PDF of 
the energy fluctuations, shown in the right panel of Fig. [1] 
(lower data set, blue circles). Again, we see a combination 
of a zero centered Gaussian and an intermittent exponen- 
tial tail. The full line is a fit to the zero centered Gaussian 
in the interval -5 < SE < 20. 

As intermittent fluctuations are rare, the temporal cor- 
relation between intermittent energy and ZFC magnetiza- 
tion changes is best observed using conditional PDF's (up- 
per curve, red squares). The conditional PDF shown only 
includes those SE values which either fall in the same or 
in the preceding small interval of width St as magnetic 
fluctuations above the threshold SM = 50. In the left 
panel of the figure, the threshold is seen to be near the 
boundary between reversible and intermittent magnetiza- 
tion changes. Note how the Gaussian parts of the full and 
of the conditional PDFs nearly coincide, while the cor- 
responding intermittent tails differ. The tail of the condi- 
tional PDF is strongly enhanced (roughly, between 10 and 
100 times) when only fluctuations near a large magnetiza- 
tion change are included in the count. A similar situation 
is observed at any low temperatures. 
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Fig. 3. (Color on line) For eight different temperatures T — 
0.2, 0.3 . . . 0.9 (circles, dots, squares, plusses, diamonds, stars, 
pentagrams, and hexagrams) the estimated cumulative prob- 
ability Prob(log tq/tq-i > x) is plotted versus x, with t q , 
q — 1,2... denoting the times at which the quakes occur. The 
quakes are identified as described in the main text using the 
filtering parameter value / = 7. The dotted lines are fits to 
the exponential exp(— ax). The collapse of the eight data sets 
demonstrates that the temperature dependence is very weak, 
except at the lowest T. In the insert, the fitted a values are 
plotted versus T. The circles pertain to the data shown in the 
main plot, and the squares pertain to corresponding data ob- 
tained using / = 8 as a filter parameter. 



Figure [2] details the dependence of the PDF of the 
magnetization fluctuations on St and T: The three PDFs 
shown in the left panel of Fig. [2] as green circles, blue 
squares and magenta diamonds are obtained in the inter- 
val [1000,3000], for T = 0.4 and for St = 20,40 and 60, 
respectively. The black line is a zero centered Gaussian fit- 
ted to the PDF in the range -80 < SM < 10. The fit is ob- 
tained by optimizing the variance of the Gaussian distribu- 
tion, Oq. The latter quantity should not be confused with 
the variance of the full distribution, which is much larger 
and heavily influenced by the tail events excluded from 
the Gaussian fit. Importantly, since the Gaussian part of 
the PDF is independent of St quasi-equilibrium ZFCM 
fluctuations at T = 0.4 can be treated as uncorrelated for 
times larger than St = 20. Secondly, since the Gaussian is 
centered at zero, any changes in the average magnetiza- 
tion are exclusively due to the intermittent events. Similar 
conclusions were reached for the spin-glass TRM magneti- 
zation, [16] and for the energy outflow in the EA model, [8] 
and in a p-spin model with no quenched randomness. [15j 
The insert describes the dependence on St of the magneti- 
zation change (SM) averaged over the observation interval 
[1000,3000]. E.g. the data point corresponding to St = 20 
depicts the average of the lowest PDF plotted in the main 
figure. The right panel of fig. [2] shows a plot of the ratio 



Oq/T versus T for the Gaussian part of the magnetic fluc- 
tuations. As the Fluctuation-Dissipation Theorem (FDT) 
applies to the equilibrium-like part of the dynamics, the 
ordinate can be interpreted as the (gedanken) linear mag- 
netic susceptibility of metastable configurations. The over- 
all shape of the curve is reminiscent of the experimental 
T dependence of the ZFCM below T g . [25J 



3.2 Temporal distribution of intermittent events 

The PDF of intermittent fluctuations in the EA model 
scales with the ratio 5t/t, where t is the age of the sys- 
tem at the beginning of the sampling interval [8j. The 
rate of quakes accordingly decays as the inverse of the 
age, in agreement with the claim that the the number 
of independent quakes in an interval (t, t') has a Pois- 
son distribution with average (nj) = a\og(t/t') The pa- 
rameter a characterizing the average is interpreted as the 
number of thcrmalizcd domains contributing in parallel to 
the fluctuation statistics. This quantity is expected to be 
temperature independent and linearly dependent on sys- 
tem size [15]. Using that the differences r q = log(t g ) — 
log(t g _i) = log(tq/t q -\) are independent random num- 
bers, all exponentially distributed with the same average 
(T q ) = 1/a, identifying the times t q , q = 1, 2, . . . at which 
the quakes occur allows one to check the temporal aspects 
of statistics (as opposed to the distributionof quake sizes) . 

The identification entails some challanges of statistical 
and patter-recognition nature. The residence times t q — 
t q -i expectedly have a broad distribution [10], which, in 
connection with a finite sampling time creates a negative 
sampling bias on large t q values |llj , even with the quakes 
unabiguously identified. Secondly, correlations may well 
be present between large and closely spaced spikes in the 
signal which collectively represent a change of attractor. 
Last but not least, identification of the t g 's within a trace 
must rely on the negative sign of the fluctuations and 
on their 'sufficiently large' size. I.e. by definition suffi- 
cient energy must be released to make the reverse process 
highly improbable within an observation span stretching 
up to 2t for a quake which occurs at t. 10J A related 
and generic property of aging systems (also directly ac- 
cessible by an intermittency analysis [5]) is that reversible 
equilibrium-like fluctuations dominate the dynamics on 
time scales shorter than the age t. A simple empirical cri- 
terion for distinguishing reversible fluctuations from irre- 
versible quakes uses an upper bound on the probability, 
P rev (t, SE) that at least one thermal energy fluctuation 
of positive sign be among the t/St observations. For small 
exp(—5E/T), this happens with probability P rev (t, SE) w 
t/5tcexp(—5E/T), where the positive number c is un- 
known. A negative energy fluctuation is labeled as a quake 
if its absolute value \SE\ satisfies P rev (t, \SE\) < 10 _n , 
where n is a positive number. For a thermal fluctuation 
which occurs at age t, this criterion produces the inequal- 
ity 



\8E\ 

quake 

f 



>/ + logt; where / = nln(10) + ln(c/<ft). 

(2) 



Paolo Sibani: Linear response in aging glassy systems, intermittency and the Poisson statistics of record fluctuations. 5 



The formal definition of the 'filter' parameter / on the 
right hand side of the inequality contains the unknown 
parameter c as well as a free parameter n. We must there- 
fore treat / as a free parameter regulating how strict a 
filter the fluctuations must pass to qualify as quakes. 

For temperatures T — 0.2, 0.3 . . . 0.9, a range spanning 
most of the EA spin-glass phase, simulations were per- 
formed in the interval t S (200, 10200). Within each tra- 
jectory, quakes are identified by applying Eq. [2] and the 
quantities r q = log(t q /t q -i) are calculated and binned. 
The procedure is repeated for 2000 independent trajec- 
tories for each value of the temperature. The magnetic 
field values H — 0.1 for t > t w = 200 and zero other- 
wise, the same as in all other simulations, were used for 
simplicity. Note however that, as also demonstrated by 
Fig. 21 the magnetic field has no discernible effects on 
the energy statistics. The standard deviation, er(i), of the 
empirical probability of a point falling in the i'th bin is 
cr ( i ) = \f{Pth{i)/N obs ), where p t h(i) is the corresponding 
(unknown) theoretical probability. N b s is the number of 
r q values collected. Replacing p t h with the corresponding 
empirical probability, yields an estimate for a(i). The em- 
pirical cumulative distribution Prob(r g = log(t q /t q -i) > 
x) is fitted to the exponential be ax , with fit parameters 
obtained by mimizing the sum of the square differences 
to the data points, each weighted by the reciprocal of the 
(estimated) variance of the data point. The prefactor b has 
no physical significance. It would be unity, where not for 
the fact that the lowest x value, for which the ordinate is 
equal to one, is located at half bin size, rather than at zero. 
To check the influence of / on the form of the distribu- 
tion, all the simulations were done twice, using the values 
/ = 7 and / = 8. The results shown in the main panel of 
Fig. [3] pertain to the / = 7 case. The results for / = 8 are 
very similar, i.e. both cases produce an exponential dis- 
tribution. The scale parameter a lacks any significant T 
dependence, except for the lowest temperature, and seems 
only weakly / dependent. From the statistics, the number 
of domains can be estimated as a « 20. Dividing the total 
number of spins by a an order of magnitude estimate is 
obtained for the number of spins in a thermalized domain. 
The corresponding linear size comes out around 6 spins, a 
figures which compares well to the range of domain sizes 
observed in domain growth studies of the EA model [23] . 
The scale invariance of the energy landscape implict in 
the T independence of a is to a large degree confirmed. 
It is expected that the largest deviations be found at the 
lowest temperatures, where the discreteness of the energy 
spectrum begins to make itself felt. 

3.3 Average energy and linear response 

The drift of the average energy /i e and average magne- 
tization fizFCM was seen to mainly depend on the the 
number of quakes nj which fall in the relevant observa- 
tion interval. For an interval (t, t') this number was shown 
to have a Poisson distribution with average 

(n/) = a\n(t/t'). (3) 



For a final check of the applicability of record dynamics we 
now compare available analytic formulas |13|16j for the av- 
erage energy, magnetic response and magnetic correlation 
function to the spin glass data. The subordination hypoth- 
esis used to derive the formulas fully neglects the effects 
of pseudo-equilibrium fluctuations. Then, all observables 
have a rather simple and generic nj dependence — a su- 
perposition of exponential functions. Averaging over the 
distribution of nj produces a superposition of power-laws 
terms, some of which can be further expanded into a log- 
arithmic and a constant term. The independent variables 
are t and t/t w for the energy and the magnetization, re- 
spectively, since all quakes contribute to the energy decay, 
while only those falling between t w and t contribute to the 
response. 

Because pseudo-equilibrium fluctuations are excluded, 
deviations between predicted and observed behavior are 
to be expected. Indeed, as discussed later, the linear re- 
sponse, data have a small additive t w dependence which 
is beyond the reach of the description and which is in 
most studies explicitely split off as a 'stationary contribu- 
tion [25] , 

The average energy and linear response are calculated 
using 4000 independent trajectories. The average energy 
fj- e (t) is found to decay toward an (apparent) asymptotic 
limit eg, according to the power-law 

Me (t) -£G = a e t x '. (4) 

In the first five panels of Fig. [4] the full line corresponds 
to Eq. 2] and the symbols are empirical estimates of the 
average energy. Each panel shows data obtained at the 
temperature indicated. For each T, the five data sets dis- 
played are for t w = 50,100,200,500 and 1000 (green, 
blue, magenta, black and red, respectively). The lack of 
a visible t w dependence shows that the magnetic contri- 
bution to the energy is negligible, as expected. The last 
panel of Fig. [4] summarizes the dependence of ec and 
of the decay exponent A e on the reduced temperature 
T/T c . The estimate T c = 0.98 [30) is used for the criti- 
cal temperature of the EA model. Interestingly, the fitted 
value of ec is nearly independent of T, and is close to 
the estimated ground state energy of the EA model, [31] 
e G = -1.7003 ± 0.008. Eq. H concurs with the observa- 
tion that a 'putative' ground state energy of complex op- 
timization problems can be guessed at early times. [32 33 
The prefactor a e has a negligible T dependence. In con- 
trast, the exponent A e has a clear linear dependence A e = 
-0.05- 0.38^. 

For A e In t < 1 , the average energy decay is logarith- 
mic, /J, e (t) ~ eg + a e + a e X e ln(t) . . ., and the decay rate 
falls off as the inverse of the age [8 15J. The coefficient 
o e A e , which receives its temperature dependence from A e , 
is proportional to the average size of a quake [15] . Even 
though quakes are strongly exothermic, their average size 
may slightly increase with T. This effect of thermal ac- 
tivation on the size of quakes is analyzed (for a different 
model) in Ref. [15] . 
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Fig. 4. (Color on line) The first five panels depict, for the temperatures indicated, scaling plots of the energy per spin, with 
the constant ec, subtracted. The values of t range from i = t w to 400t w , with t w — 50,100,200,500 and 1000, (green, blue, 
magenta, black, and red circles respectively). The black line is a fit to EqU The value of t w is irrelevant for the energy decay. 
The last panel shows the T dependence of the parameters of Eq.|4]on the reduced temperature T/T c . 
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Fig. 5. (Color on line) The first five panels display, for the temperatures indicated, scaling plots of the magnetization with a 

small t w dependent quantity b(t w ) added. The abscissa is t/t w , with t w equal to 50, f00, 200, 500 and 1000. The values of t range 

from t — t w to t — 400^. The data are fitted to the sum of a logarithm and a power-law (see Eq(5}. In the inserts, the relaxation 

rate S is plotted versus the observation time t — t w . The values of S are obtained by applying Eq. [6] to the fitted functional 

form for \izfcm given in Eq[5] The last panel shows the dependence of the parameters of Eq. [5] on the reduced temperature 
t It tu„ + „f \ ;„ „«.„«u;„„l „„„„„„„ 
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An excellent parameterization of the average ZFCM 
time dependence is given by 

HZFCAt(t/t w ) =b + a m In (j~^J + b m ' ( 5 ) 

In the theory, the left-hand side should only depend on 
t/t w , i.e. b should be constant. However, a small t w de- 
pendent term, the well known stationary part of the re- 
sponse, is present. To be able to nevertheless plot the data 
versus t/t w we shift them vertically in a t w dependent fash- 
ion, i.e. we subtract the stationary term. The magnitude 
of the shift increases with t w and with T, reaching, at its 
highest, approximately 15% of the value of 60 which fits 
the tw — 50 data. 

The first five panels of Fig. [5] show the average ZFCM 
(dots) and the corresponding fits (lines) as a function of 
t/t w , for t w = 50,100,200,500 and 1000. The simula- 
tion temperature used is indicated in each case. Sets of 
data corresponding to different t w values are color coded 
as done for the average energy. The black line is given 
by Eq. [5] with its three parameters determined by least 
square fits. The field switch at t w has a 'transient' ef- 
fect described by the power-law decay term: When t/t w 
is sufficiently large, the ZFCM increases proportionally to 
\n{t/t w ), and hence at a rate decreasing as a m /t. In the 
last panel of Fig. the decay exponent A m and the con- 
stants a m and b m are plotted versus T. The full line shows 
the linear fit — A m = 1.10 + 0.33T/T C , the dotted lines are 
guides to the eye. Importantly, the constant of propor- 
tionality a m is practically independent of T, a property 
also shared by a e . This strongly indicates that the loga- 
rithmic rate of quakes, a (see Eq. [3]) must similarly be T 
independent. 

The derivative of hzfcm with respect to the logarithm 
of the observation time, the so-called 'relaxation rate', can 
be written as 

S(t w ,t) = (t - t w ) d ^ Z ^ CM =(t- t w )r Z FCM, (6) 
at 

where Tzfcm l& the rate of magnetization increase. The 
plots of S versus the observation time t b s — t t^ are 
produced via Eqs. [5] and and dispayed in the inserts of 
the first five panels of Fig. [5] One recognizes the charac- 
teristic peak at t (, s = t w [22134) . The limiting value for 
t 3> t w is a m , precisely the asymptotic logarithmic rate of 
increase of the ZFC magnetization.. 

The sum of the TRM and ZFCM is the field cooled 
magnetization (FCM). Unlike the ZFCM, both FCM and 
TRM have a a non-linear H dependence. [34] The latter 
does not seem to affect the intermittent behavior, at least 
if numerical simulation data and experimental data can 
be treated on the same footing: Experimental TRM decay 
data [16] have the same general behavior as the ZFCM 
data just discussed. However two power-law terms, rather 
than one, are needed to describe the, more copmlex, ex- 
perimental transient. The exponents have ranges similar 
to A m , and a somewhat stronger variation with T. The 
large t/t w asymptotic behavior is a logarithmic decay, 



a\n(f/t w ), where, importantly, the coefficient of propor- 
tionality a is T independent, except very close to T g . 

With an eye to the following discussion, we recall that 
the average configuration autocorrelation function of the 
EA model, normalized to unity at t = t w , decays between 
times t w and t as 

mc (t/*™) = (*/*«) Ac ; 1 ^ *»> ( 7 ) 

where the exponent can be fitted to the linear T depen- 
dence A G (T) = —0.25^-. [13] The algebraic decay follows 
from the record dyanmics. As for the response, the the- 
oretical arguments neglect the effect of quasi-equilibrium 
fluctuations for t m t w . Furthermore, they fail near the fi- 
nal equilibration stages. Note that in the notation of Ref. [12] 
the ratio t/t w is written as 1 + t/t w . 

4 Discussion and Conclusions 

In this and a series of preceding papers |10ll2l8ll3ll5ll6j , 

we have argued that non-equilibrium events, the quakes, 
are key elements of intermittency and non-equilibrium ag- 
ing. The approach takes irreversibility into account at 
the microscopic level, stressing that thermal properties 
alone, e.g. free energies, are inadequate to explain non- 
equilibrium aging. These properties remain nevertheless 
important in the description, since the record sized (pos- 
itive) energy fluctuations which are assumed to trigger 
the quakes are drawn from an equilibrium distribution. 
As earlier discussed, [TS] the energy landscape of each do- 
main must be self-similar in order to support a scale in- 
variant statistics of record fluctuations. In this optics, the 
scale invariance of the 'local' energy landscape attached to 
each domain, rather than the scale invariance of real space 
excitations, is at the root the slow relaxation behavior of 
aging systems [35136137] . 

A different approach to non-equilibrium aging consid- 
ers a generalization of the Fluctuation Dissipation Theo- 
rem (FDT) [38 39 26 40 . Its relation to record-dynamics 
is briefly considered in below. Out of equilibrium, the FDT 
is never fulfilled exactly [H] , nor is it generally possible to 
write the linear response as a function of the conjugate au- 
tocorrelation alone. Nevertheless, for time scales t—t w << 
t w , the drift part of aging dynamics is negligible, and the 
FDT docs, in practice, apply [27]. As mentioned, conju- 
gate response and autocorrelation functions are naturally 
divided into stationary and non- stationary parts, which 
pertain to the pseudo-equilibrium and off-equilibrium ag- 
ing regimes, respectively. Adopting the standard notation 
C for the correlation and R for the magnetic response, 
the FDT reads (1 - C(t))/T = R(t)/H. For equilibrium 
data, a plot of C versus x = R/H yields a straight line 
with slope — 1/T. For aging data, the same plot produces 
a straight line in the quasi-equilibrium regime, e.g. at 
early times where C is close to one. Accordingly, a mea- 
sure of the deviation from quasi-equilibrium is the Fluc- 
tuation Dissipation ratio (FDR) [38], which is defined as 
X(C) = -Tdx(C)/dC, In the case of the EA model, the 
relation In (7 = \ c log(t/t w ), which follows from inverting 
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Eq. [JJ can be used to find x(C) an d to calculate X(C). 
The latter is nowhere a linear function of C , which is ex- 
pected, as the stationary fluctuation regime is excluded 
from the description. The effective temperature [40] is usu- 
ally defined from the large t and t w asymptotic limit of 
the Fluctuation Dissipation ratio. Effective temperatures 
may depend on the choice of conjugate observables [40] , 
and are not easily measured experimentally ,'!!) f] , but of- 
fer nevertheless a simple and appealing characterization 
of aging dynamics. 

In a record-dynamics context, fluctuation-dissipation 
like relations arise out of equilibrium because correlation 
and response are both subordinated to the same quakes. 
Due to the monotonicity of their t/t w dependencies, each 
of the two can can be written as a function of the other, 
and a FDR can thus be constructed. Asymptotically, both 
correlation and response may have an approximate log- 
arithmic dependence on t/t w : From Eq. [5l we see that 
R{t/t w ) a m ln(t/t w ) for t/t w 3> 1. Considering that 
A c vanishes for T — > 0, we also obtain, for the same 
range of t/t w values, the inequality \X c (T)\\n(t/t w ) < 1. 
When the inequality is fulfilled, Eq. [7J can be written as 
1 - C(t/t w ) « -A C (T) \n(t/t w ) cx R(t/t w ). A relation for- 
mally similar to the FDT holds, with T replaced by 

_ -\ C H _ T 
eS ~~^T~ 4T c (a m /H)- (8) 

As a m ac H, there is of course no H dependence. Im- 
portantly, T c ff cx T, ctS Ojjyi IS independent of T. For the 
parameter values obtained from the fits, T e ff > T, even 
though a general argument to support the inequality is 
lacking at the moment. 

The approach used in this paper and in ref. 16J should 
be generally applicable to check for the presence of record- 
dynamics features in intermittent fluctuation data. Clearly 
the ability to perform calorimetry experiments is essential 
to directly check the temperal statistics of the quakes. In 
the absence of such data, one may assume that quakes 
produce large intermittent changes in other observables, 
and use their statistics instead. 
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